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Electronic properties and atomic structures of W, N, S, W/N, and W/S dopings of anatase Ti02 
have been systematicaUy investigated using the density functional theory (DFT) . The exchange and 
correlation effects have been treated with Heyd, Scuseria and Ernzerhof (HSE) hybrid functional. 
Mixing traditional semi-local and non-local screened Hartree-Fock (HE) exchange energies, the HSE 
method corrects the band gap and also improves the description of anion/cation derived gap states. 
Enhanced charge carrier dynamics, observed for W/N codoped titania, can partly be explained by 
the passivative modifications of N 2p and W 5d states on its electronic structure. Following this 
trend we have found an apparent band gap narrowing of 1.03 eV for W/S codoping. This is due to 
the large dispersion of S 3p states at the valance band (VB) top extending its edge to higher energies 
and Ti-S-W hybridized states appearing at the bottom of the conduction band (CB). W/S-Ti02 
might show strong visible light response comparable to W/N codoped anatase catalysts. 

PACS numbers: 71.20.Nr, 71.55.-i, 61.72.Bb 
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I. INTRODUCTION 

In the growing area of renewable energy, Ti02 (titania) 
is the most appropriate metal oxide for photocatalytic 
processes due to its powerful oxidation and charge trans- 
port properties along with its abundance, nontoxicity, 
and stability to corrosion. It is promising in applications 
such as photogeneration of hydrogen from water, dye sen- 
sitized solar cells (DSSC), degradation of pollutants un- 
der visible light irradiation and production of hydrocar- 
bon fuels»i"— Particularly, anatase phase has received an 
increased attention for it exhibits higher catalytic activ- 
ity relative to rutile and brookitei^ Intrinsic wide band 
gap of pure Ti02 (~3.2 eV for anatase^ and ~3.0 eV 
for rutile^) confines its photon absorption to ultraviolet 
(UV) region severely limiting solar energy utilization to 
^5%. Great efforts have been made to modify the elec- 
tronic properties of titania in order to extend its optical 
absorption edge into visible region and enhance its pho- 
toresponse. For this purpose, numerous studies have pro- 
posed doping of Ti02 with substitutional cations and/or 
anions as an effective approach.—"— 

Nitrogen doped titania is considered to be one of the 
most effective photocatalysts. Although it has been ex- 
tensively studied by both experimental and theoretical 
methodsfi^"— proper description of midgap states is still 
desirable. For instance, Asahi et al,^ concluded that N 
substitution for O causes a band gap narrowing due to 
mixing of 2p states of N and O. On the contrary, Irie et 
alJ^ have suggested that the visible-light response in N- 
doped titania may be due to N 2p states isolated above 
the valence band maximum (VBM). The nature of N- 
induced modifications to the electronic band structure 
depends on the doping content i^iS^ Then, another ques- 
tion relates to increasing charge trapping rate with in- 
creasing N-doping which is associated with a localized 
state below the CBj^ Most of the theoretical studies 
failed to put forward the existence of such a trap level. 



Sulphur is another effective dopant in modifying the 
electronic structure of titania. S can be substituted ei- 
ther at an O site as an anion or at a Ti site as a cation, 
depending on the incorporation techniques.—"— Both of 
them exhibit high photocatalytic activity under visible 
light <^ Moreover, visible light absorbance increases with 
the atomic percentage of S-dopant.— i^ A reliable expla- 
nation based on the effect of S 3p states on the electronic 
structure of Ti02 is desirable. 



Monodoped Ti02 is usually inefficient for light har- 
vesting. For instance, increasing concentration levels 
of transition metal doping, such as W incorporation, 
has detrimental effects on photogenerated charge car- 
rier mobility.—"— Well localized isolated energy levels 
might act as electron-hole recombination centers reduc- 
ing photocatalytic efficiency. To passivate such trap 
states, n-p type codoping with anion and cation pairs 
has been proposed.—"— Suitably chosen codopants not 
only narrow the band gap, but also serve to counter- 
act the presence of recombination centers)^"— For this 
purpose, many experiments have recently focused on the 
additional substitution of W into N-doped titania and re- 
ported significant enhancement of catalytic activity un- 
der visible-light irradiation relative to monodoped refer- 
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We used screened exchange hybrid DFT method for a 
proper description of the electronic structure of N, S, W, 
W/N and W/S doped anatase Ti02. We attempt to an- 
alyze densities of states (DOS) in relation to visible light 
absorbance and charge carrier trapping, by comparing 
with experimental observations where explanations are 
still needed. We also discuss atomic structures, ther- 
modynamical energetics and charge states of these ions 
substituted into the anatase lattice. 



II. COMPUTATIONAL METHOD 

The spin-polarized hybrid density functional theory 
calculations were carried out using the Vienna ab-initio 
simulation package (VASP)4i Ionic cores and valence 
electrons were treated by projector-augmented waves 
(PAW) method42ii^ Plane wave basis set was used to 
expand the wavefunctions up to a kinetic energy cutoff 
value of 400 eV. We used fine FFT grids with high pre- 
cision settings throughout the calculations. 

Standard DFT usually employs local or semilocal ap- 
proximations to the exchange-correlation (XC) energy. 
This leads to erroneous descriptions for some real sys- 
tems like transition metal oxides. One way of overcoming 
this deficiency is to use hybrid functionals, where a por- 
tion of the non-local Hartree-Fock (HF) type exchange 
is admixed with a semilocal XC functional. We used the 
Heyd-Scuseria-Ernzerhof hybrid functional (HSE06)^"— 
which adopts a screened Coulomb potential. Hence, 
the exchange hole becomes delocalized around a ref- 
erence point but not beyond, reducing self-interaction 
error. Moreover, resulting rapid spatial decay of HF 
exchange improves the convergence behavior of self- 
consistent procedures. The HSE exchange is derived from 
the PBEO^^i^ exchange by range separation and then by 
elimination of counteracting long range contributions, as. 
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where a is the mixing coefficient and uj is the range sep- 
aration parameter. A consistent screening parameter of 
1X1=0.2 A~^ is used for the semilocal PBE exchange as 
well as for the screened non-local exchange as suggested 
for the HSE06 functional.— The choice of exact exchange 
contribution coefficient is important. Becked derived a 
value of 20% by fitting to atomization energy data of a 
large number of molecular species. Later, it was sug- 
gested to be 25% by Perdew et al^ We adjusted this 
mixing to be 22% with which band gaps and lattice pa- 
rameters for both anatase and rutile phases of titania 
can be obtained in good agreement with the experimen- 
tal data as given in Table HI 

We used a 108-atom supercell constructed by 3x3x1 
replication of the anatase unit cell, which ensures suffi- 
cient spatial separation between the periodic images of 
the impurities as shown in Fig. [1] Various dopings of 
Ti02 have been modeled by substitution of N (or S) at 
O and/or W at Ti sites. We also tried S@Ti configuration 
and found that this is energetically 1.246 cV/cell less fa- 
vorable than S@0. For geometry optimizations and elec- 
tronic structure calculations, the Brillouin zone was sam- 
pled with 2x2x2 mesh of special fc-points. Previously, 
the same grid was shown to be sufficient to study oxy- 
gen vacancy formation in anatase Ti02 using a smaller 
(76-atom) cellj^ Atomic positions and cell parameters 
were fully optimized until spatial components of resid- 
ual forces on each ion were below 0.015 eV/A performing 
coUinear spin polarized calculations. 




FIG. 1. Model bulk cell structure containing 108 atoms for 
W/N(S) (co)doped anatase Ti02. Substitutional dopants are 
shown for N(S) at an O site, and for W at a Ti site. 



We calculated (co)dopant formation energies using, 
Ef = i^dopod - £'puro - ?^Mw - mfiT^t^s) + nfiTi + m^O, 

where Edopcd and Epurc are the total energies of doped 
and pure supercells, while /j,w, Mn, MSj fJ-Ti and fio are 
the chemical potentials of the W, N, S, Ti and O species, 
respectively. The integer n gives the number of W cations 
and m denotes the number of N(S) anions. In thermody- 
namical equilibrium with the anatase phase, /^ti and fiQ 
must satisfy the relation /XTi02 — A*Ti + 2/xo. The amount 
of Ti and O in a growth environment influences their 
chemical potentials. High(low) values of /xq correspond 
to O-rich(-poor) conditions and can also be interpreted as 
Ti-poor(-rich) conditions from the equilibrium relation. 
Under 0-rich conditions, /zq is the half of the energy 
of an O2 molecule {E02), and fiTi is obtained through 
the condition /ixi — MTi02 ~ ^02 • Under Ti-rich condi- 
tions, /iTi is the energy of a Ti atom in its bulk unit cell 
(/^Ti'"^) ^^'^ f^o is calculated from the equilibrium restric- 
tion by /io = 5(/^Ti02 " MTi)- The chemical potentials of 
the dopants are extracted from their naturally occurring 
phases, fxw is calculated as nw = -E-wOa — §£'02 • Simi- 
larly, fis is obtained from the relation /is = ESO2 ~ -£"02 ■ 
For N, we used /in = ^^N2 • Calculated dopant formation 
energies are given in Table [Hi 

For the qualitative description of interatomic charge 
distributions, we used Bader analysis based on atom 
in molecule (AIM) theory. Local charge deple- 
tion/accumulation can be computed by integrating Bader 
volumes around atomic sites. These volumes are parti- 
tions of the real space cell delimited by local zero-flux 
surfaces of charge density gradient vector field. We cal- 
culated charge states of atomic species (see Table IIII[) 
using a grid based decomposition algorithm developed 
by Hcnkelman's groupi^ 



III. RESULTS & DISCUSSION 

For the pure Ti02, the calculated lattice parameters 
show remarkable agreement with the experimental data 
as presented in Table HI The calculated band gap val- 
ues using the bulk unit cells for both anatase and rutile 
polymorphs are largely corrected by the HSE functional. 
When we use 108-atom supercell for the anatase case, it 
slightly changes to 3.23 eV as indicated in Fig. |3^. The 
bottom of the CB is formed by highly dispersing Ti 3d. 
(see PDOS in Fig. |4^) They are not isolated but a part 
of the CB in consistency with the energy bands of pure 
anatase calculated with HSE06 by Yamamoto et al.t^ 
Calculated Badcr charges yield oxidation states of +2.84 
and —1.43 for Ti and O, respectively. These are closer to 
formal values relative to those obtained by PBE<^ 



TABLE L Computational results and experimental data (in 
parantheses) for the bulk material properties of Ti02. Cal- 
culations were performed using the HSE06 functional with a 
22% exact exchange contribution. 



Phase 


Lattice parameters 




Transition 


Band gap 




a (A) 


c(A) 






(eV) 


Rutile 
Anatase 


4.57 (4.59)^ 
3.78 (3.78)^ 


2.94 (2.95 
9.45 (9.50 




Z ^T 


2.97 (3.0)^ 
3.20 (3.2)^ 



N@0 doping : Substitution of single N atoms at O 
sites was confirmed by several reports J^"— Relaxed ge- 
ometry of N dopant in the crystal has been shown in 
Fig [2] HSE functional gave a nearest neighbor Ti-N bond 
length of 1.96 A while Ti-0 bond was found to be 1.93 A 
in its pure anatase phase. In addition, Ti-N bond along 
[001] is (2.07 A) slightly larger than that of the Ti-0 
being 1.97 A. Hence, N@0 doping has a little effect on 
the lattice structure. Consistenly, the calculated charge 
states for N and O are — 1.38e and — 1.43e, respectively, 
(see Table IIIip Hence, replacement of an O with an N is 
not expected to cause significant distortions in the lat- 
tice structure. In fact, N@0 has the lowest formation 
energy, 0.53 cV, under 0-poor conditions among all con- 
sidered in this study. (Table lll| When O is abundant in 
the environment it reasonably increases to 5.10 eV. 

For the explanation of visible light response of N sub- 
stituted Ti02, the experiments (other than the study 
Irie et al^) mainly agree on the band gap narrowing 
due to hybridization of N 2p and O 2p statesj -^^i^-^'^^i^^ 
This partial answer needs to be complemented with the 
possible presence of an impurity level below the CB as 
indicated by recombination characteristics of photogener- 
ated electron-hole pairs^l and variation of photoresponse 
under irradiation with different wavelengths ji2i22ii^ Us- 
ing X-ray photoelectron spectroscopy (XPS), Irie et ali^ 
suggested an isolated N 2p narrow band formation above 
the VB, which is attributed to be responsible for the vis- 
ible light sensitivity. Recently, deep-level optical spec- 
troscopy (DLOS) measurements of Nakano et al. have 
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FIG. 2. Relaxed structures of monodopants, N, S and W 
inside anatase Ti02 lattice. Bond lengths are in angstroms. 



located two deep levels at ^1.18 and ~2.48 eV below the 
CB.— The 2.48 eV band has been identified to behave 
as a part of the VB by mixing with the O 2p bands, 
which contributes to band gap narrowing. For the for- 
mer, the possibility of oxygen vacancies depending on the 
processes involved in preparation of the samples were ad- 
dressed as a potential origin of the impurity level closer 
to the CB<2ii2^i2i However, this issue is still unclear. For 
instance, empty N orbitals might be involved)^ 

From the theoretical side, the position and the number 
of defect states in the band gap show differences depend- 
ing on the type of DFT methods. For instance, local 
density approximation (LDA) predicts single N 2p state 
at the edge of VBM^ Generalized gradient approxima- 
tion (GGA) and GGA-I-U studies predicted it to be iso- 
lated above the VBM.^>^>iai4i j^^ ^^^ DFT-hU method, 
the supplementary on-site repulsion, U, acts only on Ti 
3d electrons and is specific to the system under consid- 
eration. In this study, we identified two N driven gap 
states (in Fig. HbHh) based on the screened Coulomb 
potential approach. One of them is a singly occupied 
bonding state showing dominant N 2p character and is 
at the edge of VBM causing a considerable band gap 
narrowing by mixing with O 2p states. Moreover, PDOS 
reveals an increasing contribution of N toward the up- 
per part of the VB. Consistently, the presence of single- 
atom N impurities as diamagnetic (singly occupied) cen- 
ters was previously observed in the electron paramagnetic 
resonance (EPR) studiesi^ The hybrid character of this 
state has been confirmed by many expcrimentsj ^^'^^" — 
This reflects a strong covalency of N centers with the 
Ti02 matrix as shown in the ground state total electronic 
density in Fig. [5^. From diffuse reflectance spectroscopy 
(DRS) measurements, Spadavecchia et al^ reported an 
"apparent" band gap of 2.88 eV at 0.4 N/Ti molar ra- 
tio, which matches with our calculations as depicted in 
Fig-Hh- Similarly, Livraghi et al. found it to be 2.84 cV 
from their EPR spectra<2^ 

Reports focusing on the charge carrier dynamics tend 
to offer the possibility of oxygen vacancy levels 0.75 - 1.18 
eV below the CB.— i^i^^ Since reduction of the charge 
separation efficiency under visible light excitation is sig- 
nificantly large, vacancies are believed to promote recom- 
bination of holes and electrons. Our HSE calculations 



predicted an antibonding energy level for one spin com- 
ponent characterized dominantly by N 2p isolated 0.63 
eV below the CB (see Fig. |4)i). In the visible region, 
it might well act as a trap level centered at the N im- 
purity site. The decrease in the quantum yields with 
increasing N concentration can also be associated to this 
empty N 2p energy level which does not conflict with 
possible existence of adjacent oxygen vacancies. Consis- 
tently, Spadavecchia et al^ opened up the possibility of 
charge transfer into empty N orbitals for the explanation 
of trapping mechanism. 

S@0 doping : Sulphur was shown to be incorporated 
into an O site of anatase lattice by oxidative annealing of 
TS2'2^'^'^ Cationic doping was also achieved by Ohno et 
al.^ Recently, Li et al. reported an alternative synthe- 
sis procedure for S substitution for 0<^ We relaxed both 
S@0 and S@Ti configurations on 108-atom supercell and 
found that the former is energetically more favorable by 
1.246 eV/cell. Optimized geometry of a S dopant replac- 
ing the lattice O is depicted in Fig. [2l In comparison to 
an oxygen, the valence electrons of the S dopant have oc- 
cupy more space, (see Fig. [SJd) Therefore, the Ti-S bond 
lengths become 2.19 and 2.31 A. while the correspond- 
ing Ti-0 bonds are 1.93 and 1.97 A in undopcd anatase. 
These local structural changes cause an increase in the 
formation energy (in Table |Tl| of an S dopant relative 
to that of N. The calculated value is 3.58 eV under O- 
poor conditions in reasonable agreement with a previous 
theoretical study<^ However, the charge state of S(0) 
is calculated with HSE as -1.17e(-1.43e) (in TablcHni 
whereas it was previously predicted by pure DFT to be 
— 0.28e and — 0.67e for S and O, respectively.— Nomi- 
nally S is expected to remain as S^~ similar to 0^~ in the 
crystal. Briefly, the HSE functional improves the spatial 
distribution of charge densities and estimates electroneg- 
ativities closer to the nominal values. 

Furthermore, screened hybrid functional better de- 
scribes the effects of the S anion on the electronic struc- 
ture of Ti02 where pure DFT fails. For instance, the 
PBE functional estimates a band gap narrowing of 0.45 
eV at an S concentration of 0.0139 (one S atom in a 72- 
atom cell) after applying a scissors correction of 1.4 eVi^ 
Our HSE calculations estimate the apparent band gap as 
2.24 eV that corresponds to an effective narrowing of 0.99 
eV for S@0 doping in agreement with UV-vis diffuse re- 



flectance spectra (DRS)<^i2^ As seen in Fig. |4}, a group 
of hybrid states localized at the edge of the VB resulting 
from O 2p, Ti 3d and dominantly S 3p coupling. This 
hybridization agrees well with the formation of S-Ti-0 
bonds as reported by Li et al. These hybrid states fa- 
cilitate a state-to-band transition at 2.24 eV (^554 nm) 
might explain the observed absorbance edge at around 
560 nm<^i^ The increase in the visible DRS intensities 
with S concentration can be attributed to increased exci- 
tation probabilities by larger DOS contribution and ex- 
tended dispersion of S-Ti-0 states. In addition, the mix- 
ing of S 3p with O 2p band states increases VB width 
by 0.16 eV relative to the VBM of pure Ti02. Similarly, 
overlap of antibonding S 3p orbitals with Ti 3d states 
decreases the CBM by 0.13 eV. These changes in the 
band edges allow band-to-band transitions starting from 
2.94 eV (^422 nm) . This might be useful in explain- 
ing the redshift of absorbance shoulder near UV in DRS 
data.i^Si^i^^ 

W@Ti doping : Experiments use sol-gel and mi- 
croemulsion techniques to incorporate W into Ti 
sites.— i^ Hence, oxygen is always present in the envi- 
ronment. In parallel, we calculated a moderate forma- 
tion energy of 3.30 eV for W@Ti doping under 0-rich 
conditions, (see Table |n| On the other hand, when Ti 
is abundant, it gets as large as 12.45 eV. Substitution of 
W at Ti in the anatase lattice results in sixfold coordina- 
tion as shown in Fig. [2] Four conjugate W-0 bonds are 
1.91 A while the bonds along [001] direction become 1.94 
A. These are slightly shortened relative to those of the 
undoped Ti02. In fact, the ionicity of W-0 bonding is 
lower then a Ti-0 bond as seen in Fig. [5]:. In parallel, the 
Bader charges around W®+ and Ti"'+ nominal species are 
found to be -1.40e (W@Ti) and -1.16e (Ti@pure), re- 
spectively. The disturbance on the lattice caused by the 
substitutional W is small and remains local. Although 
these values are lower than the formal charge states, par- 
tial inclusion of exact HE exchange improves them with 
respect to previous DFT results 4^ 

In addition, W-doping causes formation of occupied 
states 0.25 eV below the CB with W 5d and Ti 3d con- 
tributions (Fig. im), which manifest similar DOS char- 
acteristics with an oxygen vacancy state^i^ and also 
with those of Nb and Ta monodopings described by Ya- 



TABLE II. Calculated formation energies (eV) for W, S, N, 
W/N, W/S doping of anatase supercell with 18 Ti02 units. 



Ti-rich 



O-rich 



N-doped 


0.53 


5.10 


W-doped 


12.45 


3.30 


S-doped 


3.58 


8.16 


W/N-doped 


10.44 


5.87 


W/ S-doped 


15.81 


11.23 



TABLE III. Average charge states (e) of dopants and their 
adjecent Ti and O atoms from Bader analysis. 



Dopant 



N 
W 

s 

W/N 
W/S 



N 



W 



Ti(nn) 0(nn) 



-1.38 



-2.08 







+2.84 


-1.43 






+2.31 






+4.60 




-1.56 


.17 




+2.78 






+4.76 


+2.82 


-1.53 


.29 


+4.20 


+2.73 


-1.54 



mamoto et ali^ Removal of an oxygen leaves excess elec- 
trons which accumulate around Ti ions. In a similar man- 
ner, W@Ti introduces excess electrons, a portion of them 
accumulating around the W center. (Fig. [St) In fact, the 
oxidation states of W and Ti species in Table IIIII reveal 
that the W valence shell 5d^6s^ loses 4.60 of its charge 
while Ti 3d^4s^ gives 2.84 of it to form bonds with adja- 
cent oxygens. Therefore, the charge accumulation around 
the W center amounts to 1.40e being larger by 0.24e than 
that of Ti. Moreover, a very weak coupling between W 
and Ti d states can be seen in the PDOS analysis of these 
shallow states (Fig.|3]). The W 5d appears as a sharp dos 
peak for the majority spin component while Ti 3d con- 
tribution is broader involving both spins. The existence 
of impurity levels close to the CB was also discussed by 
expcrimenttal studiesj^i^ Yang et al. reported a red- 
shifted onset of optical absorption corresponding to an 
effective band gap of 2.73 eV for 3% W doping. In paral- 
lel, HSE predicts the defect states lying 2.81 eV above the 
VBM. Under visible light irradiation 19 at.% W-doped 
Ti02 absorbs up to ~445 nm (Fig. 6 in Reff2£]) that is in 
good agreement with our estimation of 442 nm. Further, 
increasing W content reduces photocatalytic activityi^ 
This can be ascribed to increased W Sd contribution on 
the defect states which can act as trap levels. 
W/N codoping : In an attempt to address the charge 
recombination bottleneck observed in N-dopcd Ti02, 
Gao et al. proposed an additional coupling to WOai^ 
Numerous studies have shown drastic enhancement of the 
photocatalytic performance of titania based on passiva- 
tion of the trapping mechanism by W-N codoping42r— 

We considered W and N in the 108-atom cell (Fig. [T]) 
with different initial configurations, and found that the 
W-N pair shown in Fig. [5t(c) is the lowest energy struc- 
ture. In agreement, Gao et alM- reported the genera- 
tion of N-W~0 linkage to be responsible for the increase 
in the visible light response. The interaction between 
W and N is much stronger with HSE, giving a bond 
length of 1.78 A which is significantly shorter than the 
GGA-fU estimation of 1.85 A4i W-N pair imposes small 
local distortions while maintaining the optimal anatase 
structural properties as pointed out by Kubacka et al^. 
For instance, T-N bond length being 2.00 A is slightly 
longer than 1.96 A of N@0 case (in Fig[2]). Supportingly, 
Kubacka et al. emphasized that W presence in anatase 
lattice is likely to influence the Ti-N bond. Moreover, 
screened hybrid functional predicts lower formation en- 
ergy for W/N codoping compared with previous theoret- 
ical studies. For example, we calculated E/ as 5.87 eV 
under 0-rich conditions. fTablellI| It was estimated to be 
^8 eV by GGA+U method using a similar formulation4i 

Electronically, both W@Ti and N@0 cells resulted in 
magnetic ground states with spin multiplicities of 0.8 and 
1.0, respectively. The energy of the W bd donor level (in 
Fig ID) is compatible with that of the N 2p acceptor level 
(in FigHJi). Hence, pairing of W with N in the lattice 
forms 5d-2p hybridization giving zero net magnetic mo- 
ment as seen in Fig HJi. This will, in turn, be effective 




FIG. 3. Optimized atomic structures of W/N(S) codopants in 
anatase Ti02. Side views (a),(b) through [010] and top views 
(c),(d) along [001] directions are presented for W/N and W/S, 
respectively. 



in passivation of electron-hole recombination at dopant 
centers. 

Analysis of the Bader charges (in Table IIIip reveals 
that N gets oxidized more by 0.7e relative to N@0 (see 
Fig.s[5K,d). This charge is partially transferred from the 
nearest neighbor W (0.16e by comparing with W@Ti) 
A large part is taken from the adjacent Ti atom such 
that its stoichiometric charge state is almost restored (by 
comparing Ti oxidation states in pure, NQO and W/N 
cases). In the presence of W-N pair, the oxidation state 
of the adjacent Ti is -^2.82 (in Table ^. Briefly, the 
charge states of N and W get closer to formal oxidation 
values. This leads to a mutual passivation of N and W 
driven defect states. Indeed, the PDOS contribution of 
W 5d in the isolated states below the CB of W@Ti disap- 
pears by mixing with N 2p gap states. Therefore, W/N 
codoping causes formation of isolated defect states 0.18 
eV below the CB characterized by Ti 3d-N 2p bonding. 
Consequently, HSE predicts a weak metallization for this 
nonstoichiometric codoping. 

When W codopant pairs with N, the position of N 
driven state at the VBM of N@0 essentially remains the 
same. The main difference, however, is the disappearance 
of N 2p acceptor level of N@0. with the presence of W 5d 
donor level. The CB edge shows similar characteristics 
with that of W@Ti. In Fig. Hi, Ti 3d-N 2p defect states 
lie 2.65 eV above the VB. This corresponds to an effective 
reduction of the band gap by 0.58 eV with respect to that 
of pure anatase. In previous experimental studies this 
narrowing was reported to be 0.55 eVM-^— However, pure 
GGA and GGA-hU studies found 0.2 and 1.11 eV gap 
narrowing values4Si^ Partial compensation of deficient 
charge by W/N pairing was also confirmed to improve 
the photocatalytic efficiency supporting our results42riS 
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FIG. 4. Total (left pane) and projected (right pane) densities of states (DOS) of pure and, W and/or N(S) substituted anatase 
Ti02, calculated with HSE06 functional. Dashed lines indicate the Fermi energies. Dotted line denote the VBM of pure Ti02. 



W/S codoping : Monodoping of both S and W has 
been experimentally achieved to be effective in increas- 
ing the optical absorption efficiencies of Ti02, suggesting 
the possibility of W/S codoping. This so that presence 
of S modifies the VB edge and brings defect states at the 
VBM by mixing with O 2p and Ti 3d band states. In ad- 
dition, shifting of the CBM and the 5d donor levels are 
induced by W substitution. Therefore, as in the case of 
W/N, coexistence of S and W species is expected to have 
a cooperative effect on band gap narrowing. Guided with 
this reasoning, we investigated the lowest energy config- 
uration of W and S substituted 108-atom anatase cell 
and found that W-S pairing as shown in Fig. ^p,d is en- 
ergetically more favorable than other possibilities. For 
instance, relaxed geometry of non-pairing W-O-S config- 
uration is 3.62 eV/cell higher in energy. W-S pair notice- 
ably distorts the lattice relative to W/N case. However, 
the effect of this distortion remains local. The coordina- 
tion of S with W and Ti increases W-S and Ti-S bond 
lengths. In fact, Ti-S bond being 2.19 A in S@0 gets 
slightly larger by 0.09 A by additional substitution of a 
W as shown in Fig. |3Ji. Since relatively large W-S (2.25 
A) and Ti-S (2.28 A) bonds do not fit in the anatase 
network, S atom relaxes out of the lattice site. This dis- 
tortion results in increased formation energies with re- 



spect to those of W/N codoping as presented in Table HIl 
which suggests synthesis under 0-rich conditions is more 
probable similar to the W/N. 

Bader analysis reveals that S codopant gets more oxi- 
dized (-1.29e) with respect to S@0 case (-1.17e). The 
charge states of the W (-1-4. 20e) and of the adjacent Ti 
(4-2. 73e) indicate slightly less electron depletion from 
around them, particularly, relative to the W/N case. 
Passivative codoping effects of the corresponding charge 
transfers can be seen in the DOS structures of the W/S 
system in Fig.s |3f and HJ. For instance, as a result of 
W 5d-S 3p hybridization, the PDOS peak of the W gets 
broadened inside the defect states relative to that of the 
W@Ti case. In fact, these states form as a result of the 
mixing of W 5d-S 3p-Ti 3d bonding states with largely 
dispersing Ti 3d CB bands. As in the W-doped and the 
W/N cases, the Fermi energy is pinned above the W 5d 
donor level similar to an oxygen vacancy situation i^^i^^ 
Hybridization of orbitals is not limited to these defect 
states at the CBM. S 3p orbitals also mix with the VB 
continuum. Moreoever, a group of largely dispersing hy- 
brid states (dominantly driven by S 3p orbitals) form as 
a part of the VB edge. The shift of the VBM into higher 
energies is similar to the S monodoping case. This results 
in an effective band gap narrowing of 1.03 eV which is 
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FIG. 5. Ground state total charge densities of (a) N-, (b) S-, (c) W-monodoped and (d) W/N-, (e) W/S-codoped anatase Ti02. 
Tiiese 3D plots iiave been cut through (001) plane slightly above the dopant species for visual convenience. 



much larger than that of the W/N codoping. The appar- 
ent band gap value becomes 2.20 eV allowing optical ab- 
sorption up to ~564 nm. Partial passivation of W 5d trap 
levels by S pairing might reduce electron- hole recombina- 
tion rate. Also, W 5d states are more delocalized relative 
to typical 3d or Ad orbitals and they strongly couple to 
the CB of Ti02 which might help passivate the carrier re- 
combination. Since impurity derived states at the VBM 
and CBM are fully occupied, interband transitions are 
much more effective than localized state-to-band transi- 
tions, due to a much larger intensity. Narrowing band 
gaps without creating isolated states are much more ef- 
fective for photocatalytic activityJ^i^S In this sense, W/S 
codoping is predicted to drastically redshift the visible 
light absorption onset involving strong band-to-band and 
also band-to-state transitions, which, in turn, enhances 
the photocalaytic efficiencies. 



IV. CONCLUSIONS 

We have systematically analyzed the atomic and elec- 
tronic structures of N, S, W, W/N and W/S doping 
in anatase Ti02 by the screened Coulomb potential ap- 
proach. Our HSE results can be useful in explaining op- 
tical absorption spectra of these systems evidenced by 
many experiments. In this way, for the N doped-Ti02, 
we identified an N 2p-0 2p hiybrid state at the VBM 
and an N 2p unoccupied trap level isolated 0.63 eV be- 
low the CB. This might serve as an alternative explana- 



tion to solve some of the controversial experiment find- 
ings. Similarly, we are able to offer an understanding 
of the increased photocatalytic efficiencies observed for 
W/N codoping. The pair mixing of N 2p orbitals with 
the isolated W 5d levels has been shown on the basis of 
their DOS structures. W-N pairing results in a band gap 
narrowing of 0.58 eV. As a trend, our HSE calculations 
might be useful in the interpretation of the passivative ef- 
fect of such a pairing of a non-metal with a d-band metal. 
Then, we examined W/S codoping in anatase. By com- 
paring modifications in the DOS structures we predict 
the origin of redshifts in the absorption light edge with 
W/N and W/S codoping systems. In fact, W/S codop- 
ing significantly reduces the optical absorption trcshold 
which allows light harvesting in the large part of the vis- 
ible range. Therefore, codoping of anatase Ti02 with 
W and S is expected to achieve highly efficienct photo- 
catalysis because the low lying photo-excitations involve 
transitions from S 3p to Ti 3d and W 5d bands through 
an energy difference of 2.20 eV and the positons of these 
states are compatible with the redox potentials of water. 
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